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ON A THREE-DIMENSIONAL FREE BOUNDARY PROBLEM MODELING 

ELECTROSTATIC MEMS 

PHILIPPE LAURENgOT AND CHRISTOPH WALKER 


Abstract. We consider the dynamics of an electrostatically actuated thin elastic plate being clamped at its 
boundary above a rigid plate. The model includes the harmonic electrostatic potential in the three-dimensional 
time-varying region between the plates along with a fourth-order semilinear pai'abolic equation for the elastic 
plate deflection which is coupled to the square of the gradient trace of the electrostatic potential on this plate. 
The strength of the coupling is tuned by a parameter A proportional to the square of the applied voltage. We 
prove that this free boundary problem is locally well-posed in time and that for small values of A solutions exist 
globally in time. We also derive the existence of a branch of asymptotically stable stationary solutions for small 
values of A and non-existence of stationary solutions for large values thereof, the latter being restricted to a 
disc-shaped plate. 


1. Introduction AND Main Results 

We focus on an idealized model for an electrostatically actuated microelectromechanical system (MEMS). 
The device is built of a thin conducting elastic plate being clamped at its boundary above a rigid conducting 
plate. A Coulomb force is induced across the device by holding the ground and the elastic plates at different 
electric potentials which results in a deflection of the elastic plate and thus in a change in geometry of the 
device, see Figure [T] An ubiquitous feature of such MEMS devices is the occurrence of the so-called “pull- 
in” instability which manifests above a critical threshold for the voltage difference in a touchdown of the 
elastic plate on the rigid ground plate. Estimating this threshold value is of utmost interest in applications 
as it determines the stable operating regime for such a MEMS device. To set up a mathematical model we 
assume that the dynamics of the device can be fully described by the deflection of the elastic plate from its 
rest position (when no voltage difference exists) and the electrostatic potential in the varying region between 
the two plates. We further assume that the elastic plate in its rest position and the fixed ground plate can 
be described by a region D in IR^. After a suitable scaling the rigid ground plate is located at z = —1 and 
the rest position of the elastic plate is at z = 0. If w = u(t,x,y) for t > 0 and (x, y) G D describes the 
vertical displacement of the elastic plate from its rest position, then u evolves in the damping dominated 
regime according to 

3fM + — (t + fl|| VmIII) Am = —A I x,i/, m(Ex,i/)) A e D , t>0, (1.1) 

with clamped boundary conditions 

u = dyU = 0, (x,i/)e3D, f>0, (1.2) 

and initial condition 

u{0,x,y) = u°{x,y), {x,y) e D . (1.3) 

Here we put 

Veip := {£dxXp,edyXp,d2Xp) , 
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Figure 1. Schematic diagram of an idealized electrostatic MEMS device. 


where e > 0 is the aspect ratio of the device, i.e. the ratio between vertical and horizontal dimensions, 
A > 0 is proportional to the square of the applied voltage difference, and xp = xp{t,x,y,z) denotes the 
dimensionless electrostatic potential. The latter satisfies a rescaled Laplace equation 

e^dlip + e^dyip + dlip = 0 , {x,y,z)^Cl{u{t)), f>0, (1.4) 

in the cylinder 

:= {{x,y,z) e D x (—l,oo) : —1 < z < u{t,x,y)} 
between the rigid ground plate at z = —1 and the deflected elastic plate. The boundary conditions for xp are 

i/;(f,x,y,z) = , {x,y,z) edn{u{t)), t>0. (1.5) 

In equation (11.11) . the fourth-order term /3 A^m with /3 > 0 reflects plate bending while the linear second- 
order term tAm with t > 0 and the non-local second-order term a|| VmH^Am with a > 0 and 

llVwlli := |VM|^d(x,y) 

account for external stretching and self-stretching forces generated by large oscillations, respectively. The 
right-hand side of ( 11.11 ) is due to the electrostatic forces exerted on the elastic plate and is tuned by the 
strength of the applied voltage difference which is accounted for by the parameter A. The boundary con¬ 
ditions (11.21 ) mean that the elastic plate is clamped. According to (ll.4l) - (ll.5b . the electrostatic potential is 
harmonic in the region Ci{u) enclosed by the two plates with value 1 on the elastic plate and value 0 on 
the ground plate. We refer the reader e.g. to ||5]|9][T2l|20l and the references therein for more details on the 
derivation of the model. 

Equations dl.lb - dl.SI ) feature a singularity which reflects the pull-in instability occurring when the elastic 
plate touches down on the ground plate. Indeed, when u reaches the value —1 somewhere, the region n(M) 
gets disconnected. Moreover, the imposed boundary conditions dl.5b imply that the vertical derivative 
dzXp{x,y, u{x,y)) blows up at the touchdown point (x, y) and, in turn, the right-hand side of dl.lb becomes 
singular. Questions regarding (non-)existence of stationary solutions and of global solutions to the evolution 
problem as well as the qualitative behavior of the latter are strongly related. Due to the intricate coupling 
of the possibly singular equation dl.lb and the free boundary problem dl.4b - dl.5b in non-smooth domains, 
answers are, however, not easy to obtain. 
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It is thus not surprising that most mathematical research to date has been dedicated to various variants of 
the so-called small gap model, a relevant approximation of (ll.ll i- (ll.5b obtained by formally setting e = 0 
therein. This approximation allows one to compute the electrostatic potential 


xp{t, x,z) 


1 + z 

1 + u{t, x,y) 


explicitly in dependence of u, the latter to be determined via 


dtu + ^A^u- (T + fl||VM||2)AM = -A . ^ '2 ' (^/3/)eD, f > 0, 

-h Mj 

subject to (11.2b and (11.3b . We note that this small gap approximation is a singular equation but no longer 
a free boundary problem. We refer to EHEIIITI and the references therein for more information on this 
case. 

The free boundary problem dl.lb - dl.Sb has been investigated in a series of papers by the authors 141 [TOl- 
EllISl, though in a simpler geometry assuming D to be a rectangle and presupposing zero variation in the 
(/-direction, see also ifT^ for the case of a non-constant permittivity profile. In this geometry the deflection 
u = u{t,x) is independent of y and ip is harmonic in a two-dimensional domain. In the present paper we 
remove this assumption and tackle for the first time the evolution problem with a three-dimensional domain 
Cl{ii{t)), assuming only a convexity property on D. More precisely, we assume in the following that 

D is a bounded and convex domain in ]R^ with a C^-smooth boundary. (1-6) 


A typical example for D is a disc. We shall see later on that condition dl.6b is used to obtain sufficiently 
smooth solutions to the elliptic problem dl.4b - dl.5b in order for the trace of Vei/i to be well-defined on 
a fact which is not clear at first glance as Ci{u{t)) is only a Lipschitz domain. Still it turns out 
that the regularity of the square of the gradient trace of this solution occurring on the right-hand side of 
(II.lb is weaker than in the case of a two-dimensional domain Ci{u{t)) studied in ll4l fT0ljT^fT5]l restricting 
us to study the fourth-order case /3 > 0 only, see Remark[3T|for more information. 

The following result shows that (ll.lb - (ll.5b is locally well-posed in general and globally well-posed for 
small a or small initial values provided that A is small as well. 


Theorem 1.1 (Local and Global Well-Posedness). Suppose ( 11.6b . Lete > 0, e [7/3,^), and consider 
an initial value uP G {D) such that u^{x,y) > —Ifor {x,y) G D and vP = dyuP = 0 on dD. Then, 
the following are true: 

(i) For each voltage value A > 0, there is a unique solution (u,ip) to (ll.lb - (ll.5b on the maximal 
interval of existence [0, Tm) in the sense that 

U G c([0,T„0,W2^^(D))nc((0,r^),w|(D))nci((0,T,„),L2(D)) (1.7) 

satisfies (frni-(fT3b together with 

u{t,x,y) > -1, {t,x,y) G [0,T,„) x D, 

and tp{t) G W|(n(M(f))) solves (ll.4b - (ll.5b in Cl{u{t)) for each t G [0, T^). 

(ii) If for each T > 0, there is k{T) G (0,1) such that 

' u{t)>-l + K{T) in D 

for t G [0, r„i) n [0, T], then the solution exists globally, that is, = oo. 

(iii) Given k G (0,1 /2), there exists := (k, e) > 0 such that Tm = oo and 

for t > 0, provided that A -|- fl|| II 2 < ffi* and 

||«0||^4e(oj<(2K)-l, u^>-l+lKinD. 
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(iv) If D is a disc in and ifi = iP{x,y) is radially symmetric with respect to (x,y) G D, then, 
for all t G [0,Tm), u = u{t,x,y) and ip = \p{t,x,y,z) are radially symmetric with respect 
to {x,y) G D. 

The global existence criterion stated in part (ii) of Theorem 11.11 involving a blow-up of some Sobolev- 
norm or occurrence of a touchdown is not yet optimal as the possible norm blow-up has rather mathematical 
than physical reasons. In the case of a two-dimensional domain Ci.{u) this condition is superfluous as shown 
in HD. Note that part (iii) of Theorem ll.ll provides uniform estimates on u and ensures, in particular, that 
u never touches down on 1 , not even in infinite time. 

Regarding existence of stationary solutions to (fOTi-(fT3]l for an arbitrary convex domain D we have: 

Theorem 1.2 (Stationary Solutions). Suppose ( II. 6 I 1 . Then, given e > 0 and k G (0,1), there are 
S := S{k,£) > 0 and an analytic function [A h-> Ux] '■ [0,^) —^ ^uch that {Ux,^\) is for 

each A G (0,(5) an asymptotically stable stationary solution to (ll.ll l- (ll.5b with Ux > —1 + k in D and 

Ta g w|(n(LfA)). 

We shall point out that while Theorem ll.2l ensures the existence of at least one stationary solution for a 
fixed, sufficiently small voltage value A, a recent result IfTSl yields a second one for (some of) these values 
in the two-dimensional case. 

When D is a disc in IR^, additional information on stationary solutions can be retrieved, in particular 
a non-existence result when a = 0. Roughly speaking, this additional information is provided by the fact 
that the operator — tA satisfies the maximum principle when restricted to radially symmetric functions 
(see Section|5]for more details). 

Theorem 1.3. Assume that D is a disc in and let e > 0. 

(i) Foranyx G (0,l)andA G (0.3(k. e)). the stationary solution (lli.Y i) to (ll.ll l- (ll.51 l constructed 
in Theorem U .2\ enioys the following properties: the function Ux is radially symmetric and Ux < 0 
in D while Ya is radially symmetric in the first two variables {x,y). 

(ii) Assume that a = 0. There are e* > 0 and a function A : (0,£*) —> (0, 00 ) such that there is no 
radially symmetric stationary solution {u,ip) to (ll.lb - (ll.51 l for A > A(£) and £ G (0, £*). 

The outline of this paper is as follows. In Section|2]we first investigate the elliptic problem (ll.4b - (ll.5b 
for the electrostatic potential in dependence on a given deflection of the elastic plate. The main result of 
this section is Proposition 12.II which implies that (ll.lb - (ll.5b can be rewritten as a semilinear equation for 
the deflection u only. The proof is rather involved and divided into several steps. Section|3]is then devoted 
to the proof of Theorem 11.11 while the proof of Theorem 1 1.21 is given in Section |4] Finally, in Section |5] 
we indicate how the proof of Theorem 11.31 can be carried out based on the corresponding proof for the 
two-dimensional case. 


2. The Electrostatic Potential Equation 
We first focus on the free boundary problem (ll.4b - (ll.5b which we transform to the cylinder 

n := D X (0,1). 

More precisely, \ct q > 2 be fixed and consider an arbitrary function v G IV^g(D) taking values in 
( — 1, 00 ), where 

{tu G 1V“(D) : u) = dytv = 0 on dD} , a > 1 + 1/p, 

{lu G W“(D) : a; = 0on3D}, a G (1/p, 1 + 1/p), 


W;(D), 


0 < a < 1/p , 


:= 




FOURTH-ORDER MEMS IN THREE DIMENSIONS 


5 


for p >2. We define 

n(z 7 ) := {{x,y,z) e D X (—l,oo) : -1 < z < v{x,y)} 


and consider the rescaled Laplace equation 

= 0, (x,i/,z) e n(t;), (2.1) 

1 “h Z 

xp^{x,y,z) = —— - -r, {x,y,z) Edn.{v) . (2.2) 

l+z;(x,y) 

Introducing the diffeomorphism 7u : n(v) ^ d given by 

Tv{x,y,z):= ' {x,y,z) Ed[v), (2.3) 

we note that its inverse is 

T~'^{x,y,y) = {x,y,{l + v{x,y))y-1) , {x,y,i])en, (2.4) 


and that the rescaled Laplace operator in ( 12.1b is transformed to the ti-dependent differential operator 


32, 


C,w := F + F afa, - 2F , a,3,w - 2F ,, 

. .._|Vu(x,y)|2 

{l + v{x,y)Y 


Av{x,y) 


{l + v{x,y)Y l + v{x,y)\ ^ 


dnW . 


The boundary value problem (I2.1b - (l2.2b is then equivalent to 


{Cv({>v){x,y,t]) = 0, {x,y,y)ed, (2.5) 

(pv{x,y,y) = y , (x,i/,y)e 8 Q, ( 2 . 6 ) 

via the transformation (pv = Observe that (I2.5b - (l2.6b is an elliptic equation with non-constant 

coefficients but in the fixed domain O. 


We now aim at studying precisely the well-posedness of (I2.5b - (l2.6b as well as the regularity of its solu¬ 
tions. To this end, we introduce for KG (0,1) and p >2 the set 

Sp(K) := |a; G Wp g(D) : ||a;||^ 2 ^p) < 1 /k and - 1 -|-k < a;(x,i/) for (x,i/) G D 

with 

Sp(K) := |a; G Wp g(D) : < 1 /k and - 1 -|-k < t(;(x,i/) for (x,i/) G D 

being its closure in ^(D). The key result of this section reads: 

Proposition 2.1. Suppose (11.6b . Let e > 0, k G (0,1), and y > 3. For each v G Sq{K) there is a unique 
solution (pv G 14 / 2 ( 0 ) to (I2.5b - (l2.6b . Furthermore there is Cj (k, e) > 0 such that 

Mvi - (pV2\\w^(^n) ^ d{K,e) 11^1 -^^2|Iw 2{D) ' e Sq{K), (2.7) 

and the mapping 

ge : Sq{K) L 2 (D), z; i—^ i-t£JZ|L |a,,^„(.,.,l)|2 
is analytic, globally Lipschitz continuous, and bounded. 

Several steps are needed to prove Proposition 12. II We begin with the analysis of the Dirichlet problem 
associated to Cn. 















6 


PHILIPPE LAURENCOT AND CHRISTOPH WALKER 


Lemma 2.2. Suppose (I1.6I ). Let e > 0, k £ (0,1), and q > 2. For each v G Sq(K) and F £ L 2 (n), there 
exists a unique solution 

O £ := {w £ W 2 ^(n) : iv = 0 on 30} 

to the boundary value problem 

— = F in Cl, <1> = 0 on 30 . (2.8) 


Proof. Since q > 2, the definition of S^(k) and Sobolev’s embedding theorem guarantee the existence of 
some constant Cg > 0 depending only on q and D such that, for v £ S^(k), 

l + v{x,y)>K, {x,y) e D, and ||t7||ci(D) < ^ ■ (2-9) 

We now claim that, due to ( 12.9b . the operator —Cv is elliptic with ellipticity constant i/(k, e) > 0 being 
independent of £ Sq{K). Indeed, for £ O, set 


flii(zt) := e^, 

a22{v) :=e^, 
assM ■■ = 
b2(v) := 


1 + q^\S7v\ 

(1 + v)^ 

( 






\\Fv, 

which allows us to write —Lv in divergence form; 


a23{^) = a32{v) := -e^ q (l^) ' 

■= 

b^iv) := -£ 2 ^ {x,y), 


-LvW = -dx (auM dxW + fll3(t’) dqw) - 3y {a22{v) dyW + fl23(^’) 3^^) 

- dq (a3i(v) dxW + a32(v) dyW + fl33(^’) dqZv) 

- bi(v)dxW — b 2 {v)dyW — b 3 {v)dyW . 

Denoting the principal of —Cy by —£„, that is, 

-Clw := - dx (aniv) dxW + ai3{v) dyw) - dy {a22{v) dyiv + a23iv) dyw) 

- dy (fl3i(u) dxtv + a 32 {v) dyW + a33(v) dyw) , 


and introducing the associated matrix 


/ 


V : = 


0 


£^dxv{x,y) 

{-YTltiFy)’' 


£^dxv{x,y) \ 

l + v{x,y)^ 

2 £^dyv{x,y) 

^ l + v{x,y)^ 

£^dyv{x,y) l + e^tf\Vv{x,y)\'^ 
l + v{x,y)^ {l + v{x,yf ) 


we observe that, for {x,y,q) £ Q, the eigenvalues of V are and 


}i± = -[t± V- 4d) 

with t and d given by 

2 , l+e^q^\Vv{x,y)\^ 

(l + i?(x,y))2 ' ■ (l + ii(x,y))2 ■ 
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By dZgl l. 


d __ 

^ ~ t ~ {1 + £^)k^ + le^CQ 


> 0 , 


which implies that —£zj is elliptic with a positive ellipticity constant depending on k and e only. Further¬ 
more, we infer from (12.9b and the definition of S^(k) that 


3 


E 

hM 


\(iiji^)\\Lo,{ci) + E < C2{K,e) 

i=l 


( 2 . 10 ) 


for all V G S^(k). It then follows from || 6 l Theorem 8.3] that, given v G S^(k) and F G £ 2 ( 0 ), the 
boundary value problem (12.8b has a unique weak solution <I> G ^(0). □ 


Next, for smoother functions v, we make use of the convexity of O to gain more regularity on the solution 
to ( 12 . 8 b . 

Lemma 2.3. Suppose (11.6b . Let e > 0 and k G (0,1). For each v G Sco{k) and F G L 2 {Cl), the weak 
solution <I> G W 2 g(n) to (12.8b belongs to W 2 g(Q) := H W 2 g(Q). 

Proof. Consider F G F 2 (n) and denote the corresponding weak solution to (12.8b by O G W 2 ^(0). The 
regularity of < 1 > and v ensure that 

G := F + bi{v) -I- b 2 {v) dy^ + b 3 {v) G 12 ( 0 ) . 

Since O is convex and (v) G Wi(n)forl < t, i < 3 , we are in a position to apply [S] Theorem 3.2.1.2] 
to conclude that there is a unique solution O G W|g(n) to the boundary value problem 

— = G in Q , O = 0 on dCl. (2.11) 

where is the principal part of the operator Cy. It also follows from 16] Theorem 8.3] that (12.11b has a 
unique weak solution in Due to (12.8b and the definition of G, the functions O and <I> are both 

weak solutions in to (12.1 lb and thus <1> = O G W|g(n). □ 


The next step is to adapt the analysis performed in the two-dimensional case in ifTSi Section 4] to derive 
an estimate on the W 2 (n)-norm of which is suitably uniform with respect to v. We begin with the 
following trace estimate. 


Lemma 2.4. Suppose (11.6b . Let p G [2,4]. There exists C 3 (p) > 0 such that 




p 

J.p(D) 


< C3(P)||W 


(3p-4)/2 

wi(n) 


IV 


(4-p)/2 
L 2 (n) ' 


w G W2^(Q) . 


Proof. Let iv G W 2 (n). Forp G (0,1) and {x,y) G D, one has 

+ p f \w{x,y,z)\^~^w{x,y,z) d,jw{x,y,z) dz . 

Jtj 

Integrating the above identity with respect to {x,y) G D gives 

+ p \zv{x,y,z)\^~'^\d„iv{x,y,z)\d{x,y,z). 

Jn 


( 2 . 12 ) 
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We next integrate with respect to t] e (0,1) and use Holder’s inequality to obtain 


p-1 
^2{p 


^ ll®ll 2 lkllLli)(n) + Pll''^llL,(p_i,(n)l|3'/«^llL2(n) 


We finally use the Gagliardo-Nirenberg inequality 

to complete the proof. □ 

Lemma 2.5. Suppose (ll.61 l. Let e > 0, k G (0,1), and q > 2. For each v G S^(k) and F G £ 2 ( 0 ), the 
weak solution <1> G W 2 g(n) to (12.8b belongs to the Flilbert space X(n) defined by 

X(n) := |tu G G W](n)| , 

and there is C^{K,t) >0 such that 

ll‘^llwi(n) + P’?^llwi(n) - ‘^4(K/e) (ll^llL2(n) + 11^llL2(n)) ■ (2.13) 

Proof. We first recall that, due to the continuous embedding of W^(D) in W^(D), there is a positive 
constant Cg > 0 depending only on q and D such that, for v G Sq{K), 


Co 


l + v{x,y)>K, {x,y)eD, and ||r||ci(D) <— 


(2.14) 


Consider F G L2{0.) and denote the corresponding weak solution to (12.8b by <I> G g(n). We begin 
with an estimate for O in W 2 (O) and first infer from (12.8b and the Divergence Theorem that 


/ F0d{x,y,q) = - ^£j,ci>d{x,y,q) 


— c 2 || p ||2 f 2 || p ||2 -I- II p ||2 

-c ll-fr 11^2(0)ll-'yllL2(n) + ll-''?llL2(n) 


— e' 


/ 

in 


d-rV 




dyV IVtlp 

- OdyO - J/tt- 


(1 + vY 


d{x,y,q) 


with 


dxV 


dyV 




Px ■■= dx<i> - y— —Py := dyd> - -—dyO, Py ■= - 

1 + v ' •' ^ 1 + v ' ' 1 + v 

Note that (12.14b ensures that 


l|V‘I>lli,(n) 


< c(k) 

Using ( 12.14b and Cauchy-Schwarz’ and Young’s inequalities, we obtain 


IP 11 ^ 
P^llL 2 (n) 


ll-PyllL 2 (n) 


ll^'7llL2(n) 


(2.15) 


f 2 ||p ||2 _|_p 2 ||p ||2 -l-llp ||2 

c IP:fllL2(n) ll''yllL2{n) + ll-'';llL2(n) 

dyV 


= [ F<Pd{x,y,rj) +e^ [ 
in in 


< ll-f|lL2(n)ll‘I’llL2(n) +£^ll‘I’llL2(n) 


-^Px + 


dyV 


-OPv 


Vti 


1 + v 

Loo{D) ^ 


d{x,y,r]) 

^’x||L 2 (n) + ll-PyllL 2 (n) 


- 2 L 


l^^llL 2 (n) + 11^1 


2 

yiiL 2 (n) 


+ c(K,e) 


U{n) 


+ ll‘I>l 


L 2 (n) 
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whence 


e^ll-P^llL{n) + £^ll-PyllL(n) + \\Pv\\l2{o) < c(K,e) 


iLzCn) 


+ 11^1 


L2(n) 


Combining (12.151 1 with the above estimate gives 


|0| 


wi{n) lll^llL(n) + ll‘I’llL2(n) 


(2.16) 


We now turn to an estimate on 3,y<I> in W 2 (n) which is established first for smooth functions v, the 
constants appearing in the estimates depending only on q and k. Indeed, assume first that, besides being in 
Sq{K), the function v also belongs to Soo{k') for some k' £ (0,1). Then O E W|g(n) by Lemma|23]and, 
setting 

^y:=dydn^, 

it follows from Lemma|Z6]below that 


/ dl<i>d^^d{x,y,q)= f \^^\^d{x,y,q), 

d^y^d{x,y,q) = \Zy\^d{x,y,q) . 

We then infer from (12.81) and (12.171) that 

- / F ^yd{x,y,q) = / dl^ Cy<^d{x,y,q) 
in in ' 

= £^llQx|li2(n) + ^^llQyllL2(n) + llQ'/llL2(n) 


(2.17) 


■I 

in 


In V 


/ 2| Vnp — (1 + v)Av 


{1 + vy 


where 


drV dyV 

Qx ■= - V ^ + 1 ;^'? ' 1 _|_ 


dy<^dy^d{x,y,q) , 


Qtj ■ — 


Owing to (12.14b we note that 


^^i/‘^llL2{n) - llQ^llL2{n) + llQyllL2(n) + IIQi/IIls 


1 + v 


(n) 


(2.18) 


Since 3^0 = dy (|3^<I>p) /2, using integration by parts to handle the last term on the right-hand side 
of the above identity leads us to 

II Qj: 11^2(0) + £^llQyllL2(n) + llQyllL2(n) 

(2rVzL_(i_p^^ \dy^{x,y,l)\^d{x,y) . (2.19) 

We now estimate successively the three terms on the right-hand side of (12.19b and begin with the first one 
which is the easiest. Indeed, it follows from (12.14b and Cauchy-Schwarz’ and Young’s inequalities that 

j^F^yd{x,y,q) < ||(l + u)f||L2(n)llQyllL2(n) < ^llQyll^in) +c('c)l|f|li2(n) ■ (2-20) 
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Next, introducing q' := q/{q — 1) G (1,2), we infer from Holder’s and Gagliardo-Nirenberg inequalities 
that 

|2 


f f2\Vv\^-{l + v)Av\ 2 ^ 

L[ -(TTtip- ) I'’’®' 


< 


2||V52||L^(D)||Vtl||r (D) + ||1 + t7||L„{D) II A^IIl„(D) 






= (ll3.,*lli,(n) + Il™„'t.|l5,,n|) l|3,*||g’„f L 

Using (12.181) and Young’s inequality we end up with 

f /2|Viip-(l + Z7)Az7\ 2 j/ n 

— 2II13,111,(0) + 2ll^^lll2(^^) ^ ^11132; II 1,(0) + ■ 

Similarly, since 2q' G (2,4), Holder’s and Young’s inequalities combined with (12.181) and Lemma lZ^ entail 
that 

/2| Vnp — (1 + zi)Az7' 


( 2 . 21 ) 


i 


{1 + vy 


|a,,d>(x,y,l)| d[x,y) 


1 r 
- ^ 


2||V52||L„(D)||Vll||mD) + ||1 +f||L»(D)l|A^llL,{D) ' A)\\l^,^o) 


<^('')ii 3 ,i>iiS'''ii 3 'i®iiw;,n;‘ 

^ ~ 1|2 I ^ liro l |2 I ^ 11 ^ ||2 


|('?+ 2 )/'? 


Inserting (I2.20b - (l2.22b in ( 12.19b leads us to 


( 2 . 22 ) 


^^llQ^cIll^in) + ^^llQyllL2(n) + llO'/llLin) - (ll^ll L{n) 


,01 


L2(n) 


from which we deduce, thanks to (12.16b and (12.18b . that (12.13b holds true. 

Since the estimate (12.13b does not depend on the regularity of v, the fact that it extends to all functions 
in Sq{K) follows by a classical approximation argument. □ 


It remains to prove the auxiliary result used in (12.17b which is recalled in the following lemma. 

Lemma 2.6. //'O G W|g(Q), then 

/ dl^dl<t'd{x,y,q) = / \dxd„<^\^ d{x,y,q), 

in ' in 

/ dyCi>d^0d{x,y,q) = / 13^3,^01^ d(x,i/, 17 ) . 

Proof. Since Q = D x (0,1) and D is a convex subset of R^, the projection pr 2 {D) of D onto the y-axis 
as well as the sections 

D[y] ;= {x G R : (x,y) G D} , y E pr 2 {D), 

are intervals. Moreover, the Fubini-Tonelli Theorem (together with Nikodym’s characterization of Sobolev 
spaces via absolutely continuous functions, see El Theorem 2.1.4]) implies that for a.a. y G pr 2 {D), 
the function •) belongs to W 2 (D[j,] x (0,1)) with <!>(•,]/,•) = 0 on 3(D[y] x (0,1)) since <1> G 
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W|g(n) C C(n). Thus, D[y] X (0,1) being a rectangle, we may apply Lemma 4.3.1.2] to conclude 


that 


Jd,.,x( 0,1) ' JD^y]X{0,l) 


/D[y]X(0,l) 

for a.a. y G pr 2 {D). Recalling that pr 2 {Ci) = pr 2 {D) is measurable and 

:= e R X R : {x,y,y) G Q} = x ( 0 , 1 ), 

the first assertion follows by integrating the above identity on pr 2 { 0 .) with respect to y and using the 
Fubini-Tonelli Theorem. The second assertion is analogous. □ 

To recover the full W|-regularity of the solution to ( 12.8b we need to have slightly smoother functions v. 

Proposition 2.7. Suppose (11.6b . Let e > 0, k G (0,1), and q > 3. For each v G S^(k) and F G £ 2 ( 0 ), 
the weak solution <I> G W 2 ig{Fl) to (12.8b belongs to g(O) and there is C 5 (k, e) > 0 such that 

ll'^llwKn) - 11^1^2(0) • (2-23) 

Proof. Consider F G L 2 {Ci) and denote the corresponding weak solution to (12.8b by O. Introducing 


}{x,y,y) ■.= 2e^q 


+ 1 - 


2 

e i] 


dxV 
1 + v 
1 + Vz;|' 

(1 + v)^ 

|Vt;|2 


{x,y) dxdy^{x,y,q) +2e^ q 


1 + v 


{x,y) dyd,^^{x,y,q) 


{x,y) d^^{x,y,q) 

Ac 


{x,y) dy^{x,y,q), 


{1 + vy 1 + c 

it follows from (12.8b that O solves 

dl^ + dl^ = } in n, 0 = 0 on dfl. 

Moreover, Lemma 1231 and the continuous embeddings of W^(D) in Wi(D) and W 2 ^(Q) in Le(n) guar¬ 
antee that / belongs to L 2 (n) with 

ll/llLin) ^ (ll^llL(n) + ITIlL(n)) • 

We then infer from ||8] Theorem 3.2.1.2] that O G (O) and inspecting the proof of IS] Theorem 3.2.1.2] 
along with |8| Theorem 3.1.3.1 & Lemma 3.2.1.1] ensures that 

||‘l>|lw|(n) ^ c(K,e)||/||^^(n) • 

We have thus shown that O G (fl) and satisfies 

ll‘^llw|(n) - c(k,£) (||‘I’|lL2(n) + ITIlizin)) • 

Finally, since the embeddings of (fl) in W 2 (H) and (D) in C^(D) are compact, we may proceed 
as in the proof of |4| Eq. (19)] to derive (12.23b . □ 

After these preliminary steps we are in a position to prove Proposition |2T| 

Proof of Proposition |2TJ Fore G Sq(K) and (x,y,q) G fl, we set 

|Vc(x,i/)p Av(x,y) 


fv{x,y,q) := Cvp = q 


l + v(x,yy l+v(x,y) 


Since W^(D) is embedded in 1A(^(D), the function /j, belongs to L 2 {Ci) with 

Wfvh^iCt) < C6(k,£) , 


(2.24) 
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and Proposition [TT] ensures that there is a unique solution e W|g(n) to 

—Cv^v = /u in n , Oj, = 0 on dCl, 

satisfying 

ll‘^''llw|(n) — ll/i'llL2(n) • (2-25) 

Setting (py{x,y,r]) = <^j,{x,y,rj) + tj for {x,y,r]) e O, the function (p^ obviously solves ( I2.5b -( l2.6l l while 
we deduce from ( 12.24b and ( 12.25b that 

ll‘?’^llw|(n) — C7 (k, e) . (2.26) 

We next define a bounded linear operator A{v) & £(W|g(n),L 2 (n))by 

^(z;)cD ;= -£„d>, o e w|g(n). 

Proposition 122] guarantees that A{v) is invertible with inverse A{v)~^ G £(£2(0), VV|g(n)) satisfying 

A{v)-'^ , ^ <C5 (k,£). (2.27) 

z:(L2(n),w2g(n)) 

We then note that 

-■^(^2)||£(w2g(n),L2(n)) < C8(K,e) |bi - i^2|Iw2(d) ^ e S,(k) , (2.28) 

which follows from the definition of and the continuity of pointwise multiplication 

Wi(D) ■ Wi(D) Wi(D) Loo(D) 

except for the terms involving d^Vi and dyVi, i = 1,2, where continuity of pointwise multiplication 

L^{D)-w^{n) -^L2(n), 

the latter being true thanks to the continuous embedding of W2 (O) in L(,{n) and the choice q > 3. Similar 
arguments also show that 

11/^1 -^21112(0) ^ C9(k,£) ll^’l -t^2|lw2(D) ' e Sq(K) . (2.29) 

Now, for V\, Vi G Sq{K), we infer from (12.27b and (12.28b that 

-^(z;2)-^||£(L2{n),w,yn)) ^ Ao{k,£) ll^^i - ^2 |Iw|(d) ' 

which, combined with (12.24b . (12.27b . ( 12.29b . and the observation that cpj, = A{v)~^fj; implies (12.7b . 

Since Wy^(D) embeds continuously in £4(0), pointwise multiplication 

W]/2(D) ■ W]/2(D) £2(0) (2.30) 

is continuous and hence, invoking lfT9l Chapter 2, Theorem 5.5] (since O = D x (0,1) is a bounded 
Lipschitz domain), the mapping 

S^(k) —^ £ 2(0), V -,1)1 (2.31) 

is bounded and globally Lipschitz continuous. Thanks to the continuity of the embedding of IV^(D) in 
£00 (D), the mapping 

S,{k) ^ Wi(D), t; ^ (2.32) 

is bounded and globally Lipschitz continuous with a Lipschitz constant depending only on k and e, and the 
Lipschitz continuity of gg stated in Proposition 12.11 follows at once from those of the mappings in (12.31b 
and (12.32b . Finally, to prove that gg is analytic, we note that Sq{K) is open in and that the mappings 

A : Sq{K) £(W|g(n),£2(n)) and [z; ^ /^] : S,(k) ^ £ 2(0) 
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are analytic. The analyticity of the inversion map ^ ^ for bounded operators implies that also the 

mapping 

Scj{K) (pv = A{v)~^fv 

is analytic, and the assertion follows as above from (12.311) and (12.321) . □ 

Let us point out that also maps S^(k:) into a (non-Hilbert) space of more regularity. 

Corollary 2.8. Suppose (11.61 ). Let e > 0, k G (0,1), and cj > 3. For p G [1,2) and a G (0,1/2) with 
a < {2 — p) /p, the mapping 

g, : S,(k) ^ , z; ^ ,<?>.(•,-,1)1' 

is analytic, globally Lipschitz continuous, and bounded. 

Proof. Given p G [1,2) and p G (0,1/2] with p. < (2 — p)/p, we may replace (12.30b in the proof of 
Pror)osition l2.1l bv the pointwise multiplication 

W]/2(D) • W]/^(D) -G IV^(D), (2.33) 

which is continuous according to |[T1 Theorem 4.1 & Remarks 4.2(d)], and deduce again from 1(19] Chapter 2, 
Theorem 5.5] that the mapping 

S(j()c) —^ Wp (D), V \d^(py{-,-,1)^ (2.34) 

is globally Lipschitz continuous. Moreover, due to the continuity of the pointwise multiplication 

w1(D)-W;(D)-g Wp^(D) 

for any a < p, see again fl] Theorem 4.1], the claimed Lipschitz continuity of follows at once from 
(12.32b and (12.34b . the proof of the analyticity of being the same as in Pror)osition l2.1l □ 


3. Proof of Theorem o 

We now turn to the proof of Theorem o Let us first note that, using the notation from the previous 
section and noticing that d^tpu {x, y, 1) = dy^u =0 for [x, y) G D due to cp^ (x, y, 1) = 1 by ( 12.6b . 

the boundary value problem ( ILlb -( IL5b can be stated equivalently as a single nonlinear equation for u of 
the form 

dtu + fA^u- (T-bfl||Vu||2)AM = \dyCpu{-,-,l)\^ (3.1) 

(1 + u) 

subject to the boundary conditions (11.2b and the initial condition (IL3b . To analyze this equation it is useful 
to write it as an abstract Cauchy problem to which semigroup theory applies. Let f > 0 be fixed such that 
7/3 < 4^ < 4 and consider G g(D) with > —1 on D. Owing to the continuous embedding of 
]N^{D) in ^1(0) and C(D), there are c > 1 and k G (0,1/2) such that 

||w||^ 2 (o) < c ^ ^ W^^{D ), 

and G S^{2k) with — 1/(2 k). Next, define the operator A G £(W|g(D), £ 2 ( 0 )) by 

Aw := fA^zv — tAw , zu G W’|g(D), (3.2) 


and recall that — A generates an exponentially decaying analytic 

|U-fAii , 11 .,-tAii , 11.,-tA 




ll£{W^|,(D)) + 11^ ll£{Wio(D)) + Ik II£{L2(D),W^|,{D)) 


< Me 


—cot 


f>0, (3.3) 
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for some a; > 0 and M > 1. Set Kq := k/ (Me) E (0, k) and introduce the mapping h defined by 

h{v) := -^geiv) +fl||Vz7||2 Az;, v G S 3 (ko) . 

According to Proposition ITTI h is well-defined on S^{kq) and there is a constant Ci{k) := Ci{KQ,e) with 
II^(^)IIl 2 (d) < Ci(k) (A-hfl||Vt;||^) ||t;||^ 2 (D), v eSsiKo), (3.4) 

and 

||;t(t;i) -;t(t;2)||L2(D) < Ci{k) [A + fl(||Vz;i||2 + ||Vz;2 ||2)] ll^l - ^2 |Iw|(D) ' Vi,V2 eSsixo) . 

(3.5) 

Consequently, we may rewrite (13. Il l as a semilinear Cauchy problem for u of the form 

dtu + Au=h{u), t>0, u{0) = . (3.6) 

Choosing i? > 2M|| Vu *^||2 we define for T > 0 the complete metric space 

Vr(K,!?) := it; e C([ 0 ,r],S 3 (Ko)) : sup ||Vz;(f)||f^ 2 {D) < 

I t6[o,r] 

endowed with the metric induced by the norm in C([0, T], S^{kq)). Arguing as in the proofs of ^ The¬ 
orem 1] and mu Proposition 3.2 (iii)] with the help of ||2] Chapter II, Theorem 5.3.1], we readily deduce 
from (13.3b . (13.4b . and (13.5b that the mapping A, given by 

A(t;)(f) := -f ^ ds, fe[0,T], u e C([0, T], S 3 (ko)) , 

defines a contraction on VjiK,^) provided that T := T(A, k, d) > 0 is sufficiently small. Thus A has 
a unique fixed point u in which is a solution to (13.6b with the regularity properties stated in 

(fO . This readily implies parts (i)-(ii) of Theorem ll.il To prove the global existence claimed in part (iii) 
of Theorem o we use the fact that a; > 0 in (13.3b and proceed as in the proof of ||4] Theorem 1] to 
establish that there is m* (k) > 0 such that A is a contraction on VjiK, (?) for each T > 0 provided that 
A -f fld < m^(K). Thus m is a global solution to (13.6b . whence Theorem ll.ll (iii). 

Finally, to prove part (iv) ofTheorem ll.il let D be a disc in Introducing then u as an arbitrary rotation 
of u with respect to (x, y) G D, the rotational invariance of (11.4b with respect to (x, y) £ D implies that u 
is again a solution to (13.6b and thus coincides with u by uniqueness. This yields Theorem ll.lh ivl. 

Remark 3.1. Besides the proof of Proposition 12.71 which requires a different approach, the range of the 
map g£ identified in Proposition \2.1\ and Corollar\ \2.8\ appears to be the main difference between the two- 
dimensional case considered in Ml? and the three-dimensional case considered herein. Since the space 
does not embed in W^[D) under the constraints a G (0,1/2) and a < (2 — p)/p with 
p G (1,2) stated in Corollarv \2.8\ we cannot identify A from the previous proof as a contraction on a 
suitable space when dealing with the second-order problem f = 0. 

4 . Proof of Theorem I1.2I 

Let K G (0,1) and note that Wl(D) embeds continuously in W^(D). Define the operator A as in (13.2b 
and note that, since A G £(lY|g(D), £ 2 ( 0 )) is invertible, the mapping 

ff:Rx (W|g(D)nS 3 (K)) ^lV|g(D), (A,i;) ^ v + \A-^g^{v) - a\\SJ v\\Ia-^A v 

is analytic with J^(0,0) = 0 and Dj,J^(0,0) = id ^4 Now, the Implicit Function Theorem ensures 

the existence of <5 = S{K,e) > 0 and a branch (L1 a)a6[o,^) -^(A/ U\) = 0 for all 

A G [0,(5). Denoting the solution to (ll.4b - (ll.5b corresponding to £1^ by ^ ^ 2 (^(^a))’ the pair 
(LTa, Ya) is thus for each A G (0,(5) a stationary solution to (ll.lb - (ll.5b . This proves the existence part of 
Theorem ll.2l 
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We next use the Principle of Linearized Stability ifTSl Theorem 9.1.1] as in ||4] Theorem 3] to obtain the 
following proposition, which completes the proof of Theorem ll.2l 

Proposition 4.1. Let A G (0,^(k,£)). There are C 0 Q,rQ,R > 0 such that for each initial value G 
^2^(0) with < Tq, the solution {u,ip) to (ll.ll i-( ll.5l ) exists globally in time and 

ll“(0 “ ^A|lw|g{D) + IIljID) — “ ^A|lw|g(D) ' f mo¬ 

using Corollarv l2.8l we can actually improve the regularity of U\ slightly. Indeed, Corollarv l2.Bl and the 
fact that U\ G (D) entail that the right-hand side of 

AU^ = -Age{U;,) + a\\VUx\\lAUx 

belongs to bVpg(D) for a G (0,1/2), a < (2 — p)/p, and p G (1,2). Now the invertibility of the 
operator A in (D), ^{D)) implies that U\ actually belongs to (D) for A G [0, ^(K,e)). 

5. Proof of Theorem I1.3I 

In this section, we restrict ourselves to the particular case where D is a disc in IR^ and assume for 
simplicity that D is the unit disc in IR^. 

To prove Theorem 1 1.3 1 (i) we argue as in the proof of Theorem ll.il (iv). Let A G (0, (5 (k, e)). Since 
( IL4b is rotationally invariant with respect to (x,y) G D, any rotation of U\ is again a solution to dl.ll ) 
in lV|g(D) n S 3 (k) and thus coincides with U\ by the uniqueness assertion of Theorem 1 1.21 The non¬ 
positivity of U\ then follows from lfT4l since ^^(Lr^^) < 0. 

Next, the proof of the non-existence statement in Theorem lL3l (ii) is a straightforward adaptation of that 
of mi] Theorem 1.7 (ii)] which is based on a nonlinear version of the eigenfunction technique. We thus 
omit the proof herein but mention that it heavily relies on the existence of a positive eigenfunction of 
the operator /3A^ — tA in Kb (D) associated to a positive eigenvalue yj. This result follows from lfT4l 
Theorem 4.7] and we emphasize that the positivity of is due to a variant of Boggio’s principle 0 and 
requires D to be a disc. Moreover, the assumptions that D is a disc and that the sought-for steady state ii is 
radially symmetric also guarantee that u is negative in D by 1(14] Theorem 1.4]. Finally, the radial symmetry 
plays again an important role in deriving suitable estimates for the auxiliary function U defined by 

—AU = M in D , U = 0 on dD . 

Indeed, U is obviously radially symmetric and its profile U defined by ZT(a/|xP -|- |yp) := hl{x,y) for 
(x,y) G D satisfies 

|arW(r)|<^, m{r)\<\, rG[0,l]. 

The previous bounds follow from the fact that — 1 < m < 0 in D and explicit integration of the ordinary 
differential equation solved by U. 
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